Variational Data Assimilation via Sparse Regularization

This paper studies the role of sparse regularization in a properly chosen basis for variational data assimilation (VDA) problems. Specifically, it focuses on data assimilation of noisy and down-sampled observations while the state variable of interest exhibits sparsity in the real or transformed domain. We show that in the presence of sparsity, the $\ell_{1}$-norm regularization produces more accurate and stable solutions than the classic data assimilation methods. To motivate further developments of the proposed methodology, assimilation experiments are conducted in the wavelet and spectral domain using the linear advection-diffusion equation

Carbon flux bias estimation employing Maximum Likelihood Ensemble Filter (MLEF)

We evaluate the capability of an ensemble based data assimilation approach, referred to as Maximum Likelihood Ensemble Filter (MLEF), to estimate biases in the CO2 photosynthesis and respiration fluxes. We employ an off-line Lagrangian Particle Dispersion Model (LPDM), which is driven by the carbon fluxes, obtained from the Simple Biosphere - Regional Atmospheric Modeling System (SiB-RAMS). The SiB-RAMS carbon fluxes are assumed to have errors in the form of multiplicative biases. Our goal is to estimate and reduce these biases and also to assign reliable posterior uncertainties to the estimated biases. Experiments of this study are performed using simulated CO2 observations, which resemble real CO2 concentrations from the Ring of Towers in northern Wisconsin. We evaluate the MLEF results with respect to the 'truth' and the Kalman Filter (KF) solution. The KF solution is considered theoretically optimal for the problem of this study, which is a linear data assimilation problem involving Gaussian errors. We also evaluate the impact of forecast error covariance localization based on a new 'distance' defined in the space of information measures. Experimental results are encouraging, indicating that the MLEF can successfully estimate carbon flux biases and their uncertainties. As expected, the estimated biases are closer to the 'true' biases in the experiments with more ensemble members and more observations. The data assimilation algorithm has a stable performance and converges smoothly to the KF solution when the ensemble size approaches the size of the model state vector (i.e., the control variable of the data assimilation problem

The Impact of Assimilating Precipitation-affected Radiance on Cloud and Precipitation in Goddard WRF-EDAS Analyses

High-frequency TMI and AMSR-E radiances, which are sensitive to precipitation over land, are assimilated into the Goddard Weather Research and Forecasting Model- Ensemble Data Assimilation System (WRF-EDAS) for a few heavy rain events over the continental US. Independent observations from surface rainfall, satellite IR brightness temperatures, as well as ground-radar reflectivity profiles are used to evaluate the impact of assimilating rain-sensitive radiances on cloud and precipitation within WRF-EDAS. The evaluations go beyond comparisons of forecast skills and domain-mean statistics, and focus on studying the cloud and precipitation features in the jointed rainradiance and rain-cloud space, with particular attentions on vertical distributions of height-dependent cloud types and collective effect of cloud hydrometers. Such a methodology is very helpful to understand limitations and sources of errors in rainaffected radiance assimilations. It is found that the assimilation of rain-sensitive radiances can reduce the mismatch between model analyses and observations by reasonably enhancing/reducing convective intensity over areas where the observation indicates precipitation, and suppressing convection over areas where the model forecast indicates rain but the observation does not. It is also noted that instead of generating sufficient low-level warmrain clouds as in observations, the model analysis tends to produce many spurious upperlevel clouds containing small amount of ice water content. This discrepancy is associated with insufficient information in ice-water-sensitive radiances to address the vertical distribution of clouds with small amount of ice water content. Such a problem will likely be mitigated when multi-channel multi-frequency radiances/reflectivity are assimilated over land along with sufficiently accurate surface emissivity information to better constrain the vertical distribution of cloud hydrometers

Evaluating Data Assimilation Algorithms

Data assimilation leads naturally to a Bayesian formulation in which the posterior probability distribution of the system state, given the observations, plays a central conceptual role. The aim of this paper is to use this Bayesian posterior probability distribution as a gold standard against which to evaluate various commonly used data assimilation algorithms. A key aspect of geophysical data assimilation is the high dimensionality and low predictability of the computational model. With this in mind, yet with the goal of allowing an explicit and accurate computation of the posterior distribution, we study the 2D Navier-Stokes equations in a periodic geometry. We compute the posterior probability distribution by state-of-the-art statistical sampling techniques. The commonly used algorithms that we evaluate against this accurate gold standard, as quantified by comparing the relative error in reproducing its moments, are 4DVAR and a variety of sequential filtering approximations based on 3DVAR and on extended and ensemble Kalman filters. The primary conclusions are that: (i) with appropriate parameter choices, approximate filters can perform well in reproducing the mean of the desired probability distribution; (ii) however they typically perform poorly when attempting to reproduce the covariance; (iii) this poor performance is compounded by the need to modify the covariance, in order to induce stability. Thus, whilst filters can be a useful tool in predicting mean behavior, they should be viewed with caution as predictors of uncertainty. These conclusions are intrinsic to the algorithms and will not change if the model complexity is increased, for example by employing a smaller viscosity, or by using a detailed NWP model

What eddy-covariance measurements tell us about prior land flux errors in CO2-flux inversion schemes

0.2 after 200 km). Separating out the plant functional types did not increase the spatial correlations, except for the deciduous broad-leaved forests. Using the statistics of the flux measurements as a proxy for the statistics of the prior flux errors was shown not to be a viable approach. A statistical model allowed us to upscale the site-level flux error statistics to the coarser spatial and temporal resolutions used in regional or global models. This approach allowed us to quantify how aggregation reduces error variances, while increasing correlations. As an example, for a typical inversion of grid point (300 km × 300 km) monthly fluxes, we found that the prior flux error follows an approximate e-folding correlation length of 500 km only, with correlations from one month to the next as large as 0.6

Toward reliable ensemble Kalman filter estimates of CO 2 fluxes

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/95036/1/jgrd18220.pd

An ensemble data assimilation system to estimate CO2 surface fluxes from atmospheric trace gas observations

We present a data assimilation system to estimate surface fluxes of CO2 and other trace gases from observations of their atmospheric abundances. The system is based on ensemble data assimilation methods under development for Numerical Weather Prediction (NWP) and is the first of its kind to be used for CO2 flux estimation. The system was developed to overcome computational limitations encountered when a large number of observations are used to estimate a large number of unknown surface fluxes. The ensemble data assimilation approach is attractive because it returns an approximation of the covariance, does not need an adjoint model or other linearization of the observation operator, and offers the possibility to optimize fluxes of chemically active trace gases (e.g., CH4, CO) in the same framework. We assess the performance of this new system in a pseudodata experiment that resembles the real problem we will apply this system to. The sensitivity of the method to the choice of several parameters such as the assimilation window size and the number of ensemble members is investigated. We conclude that the system is able to provide satisfactory flux estimates for the relatively large scales resolved by our current observing network and that the loss of information in the approximated covariances is an acceptable price to pay for the efficient computation of a large number of surface fluxes. The full potential of this data assimilation system will be used for near–real time operational estimates of North American CO2 fluxes. This will take advantage of the large amounts of atmospheric data that will be collected by NOAA-CMDL in conjunction with the implementation of the North American Carbon Program (NACP)